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Statistical analysis of the figure of merit of a two-level thermo- 
electric system: a random matrix approach 
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Abstract -Using the tools of random matrix theory we develop a statistical analysis of the trans- 
port properties of thermoelectric low-dimensional systems made of two electron reservoirs set 
at different temperatures and chemical potentials, and connected through a low-density-of-states 
two-level quantum dot that acts as a conducting chaotic cavity. Our exact treatment of the chaotic 
behavior in such devices lies on the scattering matrix formalism and yields analytical expressions 
for the joint probability distribution functions of the Seebeck coefficient and the transmission 
profile, as well as the marginal distributions, at arbitrary Fermi energy. The scattering matrices 
belong to circular ensembles which we sample to numerically compute the transmission function, 
the Seebeck coefficient, and their relationship. The exact transport coefficients probability dis- 
tributions are found to be highly non-Gaussian for small numbers of conduction modes, and the 
analytical and numerical results are in excellent agreement. The system performance is also stud- 
ied, and we find that the optimum performance is obtained for half-transparent quantum dots; 
further, this optimum may be enhanced for systems with few conduction modes. 
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Introduction. — Low-dimensional systems offer a 
wealth of technological possibilities thanks to the rich va- 
riety of artificial custom-made semiconductor-based struc- 
tures that nowadays may be routinely produced. The 
properties of these systems, including confinement geome- 
try, density of states, and band structure, may be tailored 
on demand to control the transport of heat and confined 
electrical charges as well as their coupling. This is crucial 
to further improve the performance and increase the range 
of operation of current thermoelectric devices. 

From a thermodynamic viewpoint, thermoelectric sys- 
tems connected to two thermal baths at temperatures T^ ot 
and T co id, use their conduction electrons as a working fluid 
to directly convert a heat flux into electrical power and 
vice- versa with efficiency r\. As for all heat engines, one 
seeks to increase 77, the maximum of which, ^ max , is related 
to the figure of merit of the system ZT 1,2: 



Vc 



VI + ZT - 1 



yi + ZT + TcoM/Ij 



(1) 



hot 



where 77c = 1 — ^coid/^hot is the Carnot efficiency. In the 
frame of linear response theory, the figure of merit may be 
expressed as: ZT = aS 2 T '/ '/€, where S is the thermopower, 
T is the average temperature across the device, and a and 
k are the electrical and thermal conductivities respectively, 
with k = n e + ftiat , accounting for both electron and lattice 
thermal conductivities. Equation (HI) thus clearly shows 
that for given working conditions, ZT is as good a device 
performance measure as the efficiency 77 is. 

Quantum dots may be defined as confined electronic sys- 
tems whose characteristic size is much smaller than the 
electron mean free path; the size and shape of the re- 
gion where electrons are confined are controlled by gates. 
Amongst the various low-dimensional thermoelectric sys- 
tems that have been studied, quantum dots keep attract- 
ing much interest because of their narrow, Dirac-like, elec- 
tron transport distribution functions [3] or, equivalently, 
sharply peaked energy-dependent transmission profiles T, 
which permit obtainment of extremely high values of elec- 
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tronic contribution to ZT 4 . There is a simple relation- 
ship between the electrical conductivity a appearing in the 
definition of the figure of merit and the transmission func- 
tion T: (J — (2e 2 /h)T, the proportionality factor being 
the quantum of conductance (e is the electron charge, and 
/i, Planck's constant). 

The treatment of confined systems may be achieved, 
without loss of pertinence, by means of very simple mod- 
els such as, e.g., the resonant level model of thermal effects 
in a quantum point contact p\ and the two- level model of a 
mesoscopic superconductor J6] , which capture the essential 
physics of the actual systems. The use of non-interacting 
two-level model systems for quantum dots or very small 
electronic cavities may be justified if these are strongly 
coupled to the reservoirs, i.e. if the confinement yields a 
mean level spacing that is large 7 compared to the charg- 
ing energy e 2 /C, C being the capacitance of the dot [8J[q] . 
These models also are justified for probing the spectrum 
edges of mesoscopic systems 10, 11 where the physics is 



different from that of the bulk spectrum 10, 12 



If one of the gates generates a random potential, the 
quantum dot behaves as a chaotic cavity; measurements of 



thermopower and analysis of its fluctuations 13 based on 
random matrix theory of transport [l4] demonstrated the 
non-Gaussian character of the distribution of thermopower 
fluctuations. In this letter, we analyze the statistics of ZT 
in a two-level system connected to two electronic reservoirs 
set to two slightly different temperatures so that the tem- 
perature difference, AT, is small enough to remain in the 
linear response regime: the voltage induced by thermo- 
electric effect, AV, is given by S = -AV/AT. The mu- 
tual dependence of the transport coefficients that define 
ZT (e.g., the Wiedmann- Franz law for metals) raises the 
problem of finding which configuration of the mesoscopic 
system may yield the largest values of ZT, and hence offer 
optimum performance. 

The Letter is organized as follows: In the next section, 
we introduce the model system, the definitions and no- 
tations we use throughout the paper, and the theoretical 
framework of our analysis. In the third section, we analyze 
the numerical statistical results obtained for the probabil- 
ity distributions of the figure of merit, thermopower, and 
power factor. We end the Letter with a discussion and 
concluding remarks. 

Model. — The two- level model presented in this Let- 
ter, is the minimal model pertinent for the description of a 
cavity with two conducting modes presenting a completely 
chaotic behavior 11 . The main assumption on the ther- 



mal properties of the system is that n is constant; this is 
justified on the condition that n e <C ftiat ^J] (^lat is ne- 
glected in the present work: our interest is in the electronic 
transport properties only). 

Basics. The two-level system tight-binding Hamil- 
tonian reads as the sum of three contributions: H = 
Hi+H s +H c where Hi = J2k e k c k c k is the Hamiltonian of 
the two (i.e., left and right) leads; c\ and ct are the second- 
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Fig. 1: (Colour on-line) Schematic figure of the two level system 
consisting of negatively charged top gates (yellow) on a two 
dimensional electron gas. The potential on the right three top 
gates may be varied randomly in order to slightly modify the 
shape of the cavity and therefore obtain a statistical ensemble 
[13] . The red (blue) part represents the hot (cold) electronic 
reservoir connected to the cavity. The inset shows the system 
as modeled by the Hamiltonian. 



quantized electron creation and annihilation operators in 
the state k. The two energy levels and their coupling are 
characterized by the Hamiltonian H s '- 



n s 



Vi V X 2 
l>21 V 2 



(2) 



where the on-site potentials vi, V2 and the coupling v 12 
are chosen to be random to describe the chaotic behavior 
of the system. We assume that the left (right) lead is 
only connected to the site v\ (^2). The contribution H c 
involves the coupling matrix W entering the definition of 
the scattering matrix [18] : 



S = l-2mW ] - 



1 



H s 



:W 



(3) 



where E is the self-energy of the two leads. To keep 
our analysis on a general level, we only give its form: 
E = A - Z7TWW 1 ", with A being the real part. With the 
assumption of symmetric reservoirs, the self-energy is pro- 
portional to the identity, in which case W = yT/27rl2x2, 
where T — — 2^E characterizes the broadening. Through- 
out the Letter, we adopt the same notations for scalars 
and their corresponding matrix form when this latter is 
proportional to the identity. 

The scattering matrix S is unitary. In the absence of 
a magnetic field, time reversal symmetry is preserved and 
the matrix S is therefore symmetric, which implies v\i = 
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i>2i. Moreover, for simplicity, we consider in the first part 
of this Letter, systems with left-right spatial symmetry 
(v\ = i;2 ), so that the reflection from left is the same as 
that from right. The matrix S may thus read: 



S 



r t 
t r 



(4) 



where r and t are the reflection and transmission ampli- 
tudes respectively. With these definitions, the transmis- 
sion of the system and the Seebeck coefficient read: 



T=\t\ 



and S 



dln(T) 
de 



(5) 



The Seebeck coefficient is obtained at low temperatures 
with the Cutler-Mott formula n9\. Here, it is expressed in 
units of |^/cgT (ommiting the sign). The transmission of 
the system can then be directly obtained using the Fisher- 



Lee formula 20 : 



t = tg 12 tg\ 



(6) 



121^12 

where the off-diagonal element of the Green's matrix is 

G12 — \ (e _ e _^"^_ e _ S n , with ei and 62, eigenvalues of 
H s . And, we may thus write: 



T(e) = 



(ei-e 2 ) S 



4 |( e - ei -S)( e -e 2 -S) 



(7) 



which is valid if the left /right symmetry is satisfied. Now, 
combination of Eqs. (pi and u?\ yields an expression of the 
Seebeck coefficient containing the scattering matrix <S: 



= a Tr / o o_ 

2 V 2i 



with a = 4|1 - E|/r, 6 = Arg(l - £), and £ = <9 e E. 
For an energy e = eo, which corresponds in general to 
the middle of the conduction band of a semi-infinite lead 
(or half- filling limit), the imaginary part of the self-energy 
derivative vanishes: 9E(eo) = 0. We stress that Eq. (|8| 
may also be derived if the constraint on left-right symme- 
try is relaxed, and that it constitutes an important result 
upon which all the subsequent statistical analysis is based. 

Statistics of the Seebeck coefficient. In this work, the 
scattering matrix S is a random variable, which we assume 
to be uniformly distributed; hence S belongs to circular 
orthogonal ensembles (COE) 21 . The statistics of the 
Seebeck coefficient thus follows: 

3 iG S-e~ iG S^ 



V e (S) 



H 



s 



Tr 



2% 



5S (9) 



The Haar measure 5S is invariant under the transforma- 
tion S — »■ VSV T for any arbitrary unitary matrix V (the 
superscript T denotes the transpose of a matrix). With 
V = e ie / 2 l 2X 2, we obtain: 



™=/W5-H^ 



5S = =P e=eo (S/a) 

a 

(10) 



where a = a/ao with ao = a(eo). Equation (10) shows 
that the general probability distribution of the Seebeck co- 
efficient at arbitrary energy can be deduced directly from 
the simpler one obtained at eo- 

Note that for ease of notations, we retain ?asa generic 
notation for all the probability distributions in the subse- 
quent parts of the letter. 

Joint probability distribution function. A similar anal- 
ysis of the transmission T shows that its distribution does 
not depend on energy 22 . Since the transport coefficient 
T and S are tightly related by the transport equations, we 
need the joint probability distribution (j.p.d.f) V e (T, S) to 
study any observable depending on both T and S. The 
starting point for this is the rewriting of the scattering 
matrix using the following decomposition: 



S = R 2 











R 



where the rotation matrix R is denned as: 
R 



V2\ 1 



(11) 



(12) 



The eigenphases 0\ and #2 are independent and uniformly 
distributed in [— 7r, +7r] which is the consequence of S be- 
ing taken from COE with the left-right symmetry 24 . 
Now, we analyze the j.p.d.f. of the Seebeck coefficient and 
transmission at the half-filling limit e = cq 25 , with: 



T = sin 2 [(0i - 2 )/2] 



(8) and 



5 = a o sin[(0i+0 2 )/2]cos[(0i 
The j.p.d.f. may now take the form: 

V e=eo {S, T) = {5(T- sin 2 u) 5 (S ■ 



0/2] 



(13) 
(14) 



ao sinvcosii)) n v 
(15) 
where two independent and uniformly distributed vari- 
ables u = 9l ~ 92 and v = ei +° 2 are introduced. Integration 
over these variables and use of Eq. ( 10 ) yield: 



Ve(S,T) = 



1 



7iVT(l - T) 7r N /a 2 (l - T) - S 2 



(16) 



which constitutes one of the main results of this work since 
it shows the possible values of the couple (5, T): V e (S, T) 
is non-zero only if the following condition is satisfied: 



a 2 {l-T)-S 2 >0, T < 1 



(17) 



The relationship between the Seebeck coefficient and the 
transmission function is thus constrained by a parabolic 
law. This may be checked by sampling the S matrices 
belonging to the COE with the left/right symmetry, from 
which the transport coefficients T and S are numerically 
computed with Eqs. (pi) and (pi). The parabolic law is 
clearly shown on Fig. (|2| where the j.p.d.f. V e (S,T) is 
plotted against T and S. 
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Fig. 2: (Colour on-line) Probability densities: l)Top left: V(T) 
for the transmission. 2) Bottom left: V(S) for the Seebeck 
coefficient. 3) Top right: j.p.d.f. of the couple (S,T) at the 
half filling limit. 4) Bottom right: j.p.d.f. of the couple (T, ZT) 
at the half filling limit. 



To study the system performance, we compute the fig- 
ure of merit ZT whose maximum is reached when the 
power factor p = TS 2 is reached. Satisfaction of this latter 
condition lies on the existence of the probability V e (T, S): 



a(l - T) - S 2 > 



olT(\ -T)> TS 2 
max(aT(l - T)) > P (18) 
a/4 > p 



We thus find that ZT reaches its maximum if the trans- 
mission T = \ . Physically, this means that the best pefor- 
mance results from the conditions that make the system 
half-transparent; this is numerically checked on Fig. ¥A. 

Marginal distributions. The marginal distributions 
V e (T) and V e (S) are obtained from the integration of 
V e (S, T) over S and T respectively: 



P e (T) 
V C (S) 



1 1 

* y/T(l ~ T) 

2 



-X 



1 



S 2 /a 2 ) 



(19) 
(20) 



where /C is the complete elliptic integral of the first kind. 
It is interesting to note that V e (T) does not depend on the 
energy since it derives from the sole assumption that the 
scattering matrix is uniformly distributed (Dyson's CE) 
with the left/right spatial symmetry [24]. As for all the 
physical observables with no energy derivative in their ex- 
pression (shot noise for example), their distributions do 
not depend on the number of levels N in the system: for 
all N > 2, we obtain the same distribution, as it can be 
seen with the decimation method [TT]|T5] as on the condi- 
tion that S G {CE}. 



The case of the Seebeck coefficient is different: its distri- 
bution V e (S) depends on the number of levels in the sys- 
tem since its expression contains an energy derivative. The 
only situation where the Seebeck coefficient, scaled with 
the appropriate parameter (a in our case), is independent 
on the number of levels is at the edge of the Hamiltonian 
spectrum [To] where the result can be obtained by assum- 
ing the simpler two-level system considered as the minimal 
cavity for the two-mode scattering problem \T\\ . 

All the results presented so far in this Letter were ob- 
tained and discussed assuming a two-level system with 
left /right spatial symmetry and S G {CE}. We must now 
see whether these hold when the assumption of left/right 
spatial symmetry is relaxed. 

Relaxation of the spatial symmetry assumption. As- 
suming that there is no left /right symmetry, the two 
eigenphases 6\ and #2 of Eq. ( pTj ) are no longer indepen- 
dent, and we must consider a distribution of the form: 
V(0 U 2 ) oc \e iQl - e i92 \ [2l]. The treatment is the same 
as for the previous case, and while we obtain a different 
j.p.d.f. for the Seebeck and the transmission coefficients: 



Ve(T,S) = 



an 2 y/T(l - T) 



IC^l-S 2 /a 2 (l-T)) 



(21) 

the mathematical constraint on the couples (T, S) to ob- 
tain a non-vanishing joint probability distribution is ex- 
actly the same as Eq. (18); this implies that the best sys- 



tem performance is obtained when it is half-transparent: 
The marginal distributions of S and T are also different: 



-Pe(S) 



T>*{T) 



In 



2VT 



\S\/a 



an \ 1 + a/1 " S 2 /a 2 



(22) 



(23) 



Equation (22) is consistent with results of Refs. 26,27 



obtained with different methods and Eq. 
the result of 



(23) confirms 



10 . We see in both distributions, Eqs. (16) 



and (21), that the variables (T,S) are not independent 



but it is interesting to note that if we define a new variable 
X = S/ y/{l — T) then we have two independent variables, 
and the j.p.d.f. becomes multiplicatively separable: 



r e (T,X) = V e (T)xV e (X) 



(24) 



where we have V e (X) = 1/iry/a 2 — X 2 in the symmetric 
case and V e (X) = -^/C(yl — X 2 /a 2 ) when the symme- 
try is relaxed; V e (T), given above, is the corresponding 
marginal distribution. It is worth mentionning that S and 
X also constitute a couple of independent variables. 

Discussion. — 

Seebeck coefficient and density of states. The local 
density of states in the central system, when connected to 
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Fig. 3: (Colour on-line) Probability density of the power factor 
p for different values of the number of conduction modes: blue 
curve M = 4, purple curve M = 8, and brown curve M = 
16. The number of energy levels in the cavity equals the total 
number of modes. In the inset, the standard deviation of the 
power factor is shown as a function of M. The results are fitted 
with the function \/5p 2 ~ 0.08/M. The numerical results are 
obtained by sampling 2.1 x 10 5 scattering matrices. 



the leads 28 , reads: 



P(e) 



-3Tr- 

7r e 



u s 



which combined with Eq. (25) yields: 



5p(e) _l m fS + S^ 
W)' 



Tr 



(25) 



(26) 



where 5p = p — p and p(e) = 2/71T. The analysis of this 
expression shows that the relative change in the density 
of states differs from the Seebeck coefficient, Eq. (pi), but 
the interesting result is that both Sp(e)/p(e) and S/a have 
exactly the same distribution, which may be seen by us- 
ing the invariance of the Haar measure under the trans- 
formation S — >• — iS. We finally add to this set of vari- 
ables the scaled Wigner time Sr w /a which is related to 
the time spent by a wavepacket in the scattering region: 
r w = -ihTi(S^dS/de) [IT]. 

Generalisation to M modes of conduction. Here, we 
generalise and investigate the performance of a system 
with a large number of conduction modes. We concentrate 
on systems with 27V levels connected to 2M independent 
and equivalent leads (M from the left side and M from 
the right side) . We suppose that the number of modes M 
is the same as the number of levels: M = N. This case is 
pertinent because it corresponds to the model of minimal 
chaotic cavities, which gives the same results as for the 
general case N > M at the edge of the spectrum 111]. We 
also assume no spatial symmetry. 

The form of the scattering matrix which is now of size 
2M x 2M remains unchanged [see Eq. pi], and the trans- 
mission is given by T = Tr(ttf). For large numbers of 
modes, the shape of p.d.f. of the transmission tends to 




Fig. 4: (Colour on-line) Standard deviation of the Seebeck co- 
efficient as a function of the number of modes for a minimal 
chaotic cavity (M = N). The Seebeck coefficient scales as 1/M 
for large numbers of modes. 



that of a Gaussian distribution for the typical values of T 
around its mean which for M large is given by T ~ 4f . 
The variance of this Gaussian is equal to (ST 2 )=l/8 which 
does not depend on M because of the universal charac- 
ter of the conductance fluctuations 14 . It is worth men- 
tioning that the tail of the distribution, which describes 
the atypical values of the conductance has a non- Gaussian 
form 30 . At first glance, this scaling of the mean trans- 
mission seems to favour an enhancement of the power fac- 
tor and hence of the figure of merit; however, our study of 
the Seebeck coefficient yields a different conclusion as we 
now explain. 

The Seebeck coefficient may be expressed by writing 
the derivative of the scattering matrix and using the ex- 
pression S = 2$l(Tr(ffl)/T) (where the dot refers to 
the energy derivative and 5ft is the real part); then it is 
computed numerically assuming modes with a self-energy 
E = 6/2 - iy/1 - (e/2)2 



31 



Still considering the half- 
filling limit since the generalisation to arbitrary energy 
presents no particular difficulty, we find that the mean 
value of the Seebeck coefficient is always zero: (S) = 
because the probability to obtain a positive thermopower 
is the same as that for a negative coefficient. The most in- 
teresting result is the standard deviation which appears to 
scale as: (SS 2 ) 1 ^ 2 = 0(1/M). This clearly demonstrates 
that increasing the number of conduction modes in the 
system yields a decrease of the thermopower as shown in 
Fig. [4] where the standard deviation of thermopower is 
represented as a function of the number of modes on each 
side (left and right). The numerical results were obtained 
by sampling scattering matrices of size 2M x 2M from the 
circular orthogonal ensemble and computing a histogram 
from which the standard deviation was obtained. 

Since p oc S 2 , the lowering of the system performance 
induced by a decrease of the thermopower, which is faster 
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Fig. 5: Gaussian fit of the distributions V(T — T) and 
V(MS/a) (left panels), and j.p.d.f. for the couples (T,ZT) 
and (T,MS/a). 



than the increase due to conductance, is a reflection of the 
fact that the typical power factor is thus estimated to scale 
as 1/M. We deduce from this analysis that the distribu- 
tions V(T — T) and V(MS/a) do not depend on M for 
large values of M. In Fig. [5] we see a very good agreement 
between the distribution V(MS/a) and the Gaussian cen- 
tered on zero. We also see from the j.p.d.f of ZT and T 
that high values of the figure of merit are obtained for 
T ~ M/2, which generalises the result obtained for the 
two-level system. 

Concluding remarks. — We investigated the quan- 
tum thermoelectric transport of nanosystems made of two 
electron reservoirs connected through a low-density-of- 
states two-level chaotic quantum dot, using a statistical 
approach. Analytical expressions based on an exact treat- 
ment of the chaotic behavior in such devices were obtained 
and analyzed to provide the conditions for optimum per- 
formance of the system. We find that the optimum effi- 
ciency is obtained for half-transparent dots; this optimum 
may be enhanced for systems with few conducting modes, 
for which the exact transport coefficients probability dis- 
tributions are found to be highly non-Gaussian. 
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